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ABSTRACT 



We estimate the rate at which the proposed space gravitational-wave interferometer LISA 
could detect intermediate-mass black-hole binaries, that is, binaries containing a black hole in 
the mass range 10 - 100 M© orbiting a black hole in the mass range 100 - 1000 M©. For one- 
year integrations leading up to the innermost stable orbit, and a signal-to-noise ratio of 10, we 
estimate a detection rate of only 1 per million years for 10M©/100M© binaries. The estimate uses 
the method of parameter estimation via matched filtering, incorporates a noise curve for LISA 
established by the LISA Science Team that is available online, and employs an IMBH formation 
rate model used by Miller (2002). We find that the detectable distance is relatively insensitive to 
LISA arm lengths or acceleration noise, but is roughly inversely proportional to LISA position 
errors, and varies roughly as T 1 / 2 , where T is the integration time in years. We also show that, 
while all IMBH systems in this mass range may be detected in the Virgo cluster up to 40 years 
before merger, none can be detected there earlier than 400 years before merger. An extended 
LISA mission that enabled 10-year integrations could detect IMBH systems at the Virgo cluster 
1000 years before merger, and systems in galactic globular clusters a million years before merger. 
We compare and contrast these estimates with earlier estimates by Miller (2002). 

Subject headings: gravitation: gravitational radiation, black holes 



In recent years, intermediate-mass black holes 
(IMBH) , holes with masses between hundreds and 
thousands of solar masses, have attracted con- 
siderable interest, both as relics of the evolution 
and structure of globular clusters, and as possible 
sources of gravitational radiation for both ground- 
and space-based laser interferometers (Miller & 
Colbert 2004). Miller (2002) has estimated the 
rate of detectable IMBH binaries by the proposed 
space antenna LISA to be given by 



1. Introduction and summary 




where h is the Hubble parameter in units of 
100 kms -1 Mpc~ , /tot is the total fraction of 
globular clusters that have IMBH, /i and M are 
the reduced mass and total mass of the binary re- 
spectively, and Af m j n and M max denote the range 
over which such IMBH may exist in globular clus- 
ters. 
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A key ingredient in calculating this rate is the 
distance to which LISA could detect binary IMBH 
inspirals for a given signal-to- noise ratio (SNR). 
Miller (2002) used estimates of signal-to-noise ra- 
tios for binary inspiral derived by Flanagan & 
Hughes (1998) from semi-analytic LISA noise 
curves. We have recalculated this distance us- 
ing standard methods based on matched-filtering 




Reduced Mass Parameter r\ 

Fig. 1. — Distances in Mpc reached by LISA for 
IMBH binary inspiral, for SNR=10 and one year 
of integration prior to merger. 

for binary inspiral, and using an up-to-date LISA 
noise curve available online. 2 We find significantly 
smaller distances and rates reached for a given 
SNR and with a significantly different dependence 
on mass than those obtained by Miller (2002). 
Figure I shows distances reached for SNR=10 
and one year of integration leading up to merger, 
as a function of total mass M and reduced mass 
parameter 77 = ■rnxm^jM' 1 (0 < 77 < 1/4). For 
M = IOOMq and reduced mass p = rjM = 1OM , 
we find a distance Dl ~ 18 Mpc, 11 times smaller 
than Miller's equation (19) (we assume that all 
masses are suitably redshiftcd). The rate of de- 
tection we find is given by 

( M max \ 13/4 / |n M m „ r' ! 

X ^ 1OOM ) \ U M min J yr ' U 

We also analyse LISA's reach for IMBH sys- 
tems at epochs earlier than the year leading up 
to merger. Figure 2 shows distances MS (X function 
of years before merger for various IMBH systems, 
for SNR=10 and one year of integration. While 
all systems in this mass range are detectable in 
the Virgo cluster (at 19 Mpc) within 40 years 
of merger, none is detectable there earlier than 

2 The Sensitivity Curve Generator was originally writ- 
ten by Shane Larson and may be found online at 
http:/ /www. srl. caltcch.edu/" shane/scnsitivity/MakeCurve 
.html 
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Fig. 2. — Distances in Mpc reached by LISA for 
IMBH binary inspiral, for SNR=10 and one year 
of integration vs. years before merger, for various 
black-hole pairs (in M Q ). 

400 years before merger, and at 1000 years be- 
fore merger, only systems closer than 7 Mpc are 
detectable. On the other hand, we find that, for 
systems in this mass range, the distance reached 
varies roughly as the square-root of the integra- 
tion time, so that an extended LISA mission that 
enabled 10-year integrations, could reach roughly 
three times farther, 

The primary reason for this downward revision 
in the reach of LISA for IMBH is our use of a more 
current sensivity curve for LISA, whose noise level 
is substantially higher than the older Flanagan- 
Hughes curve used by Miller. The rest of this pa- 
per gives the details to support this conclusion. In 
§2 we determine the distance reached by LISA us- 
ing standard SNR calculations, and in §3 we use 
the method of Miller (2002) to estimate the detec- 
tion rate of systems in the final year of inspiral. 
In §4 we discuss the results. 

2. Calculation of Distance Reached by 
LISA 

The signal-to-noise ratio p for a gravitational- 
wave signal whose Fourier transform is h(f) in a 
detector of noise spectral density S n (f) is given by 

where /j and ff are the initial and final frequen- 
cies between which the signal is integrated (see 
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Cutler et al. (1993); Cutler & Flanagan (1994); 
Finn (1992); Finn & Chernoff (1993); Poisson & 
Will (1995) for details of the method of matched 
filtering). In the "restricted post-Newtonian" ap- 
proximation, where the wave amplitude is given 
by the quadrupole approximation and the wave 
phase is given by a high-order post-Newtonian ex- 
pression, h(f) is given by 
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(4) 



where / max is the largest frequency for which the 
wave can be described by the restricted post- 
Newtonian approximation, often taken to corre- 
spond to waves emitted at the innermost stable 
orbit (ISCO) before the bodies plunge toward each 
other and merge. A useful approximation for this 
frequency is (we use units in which G = c = 1) 



/max = (6 3/2 7tMY 



(5) 



where M is the total mass of the system. After 
averaging over all angles, the amplitude A is given 

by 

1 M 5/6 
A = V30tt2/3 D L ' (6) 

where Dl is the luminosity distance of the source, 
and M. — if^M is the "chirp" mass. In general 
relativity, the phasing function x &(/) can be given 
to high post-Newtonian order, but is not relevant 
for our purposes. Combining equations (3), (4) 
and (6), we obtain the expression for Dl 



D, 



2 M 5/e 1 



15 p 



r 2/3 



B{7) 



1/2 



where 



B(7) = 



ff f-7/3 



fi 



Sn(.f) 



df. 



(7) 
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For quasi-circular binary inspiral, the initial and 
final frequencies can be related to the integra- 
tion time T using the expression obtained from 
gravitational radiation damping of the orbit at 
quadrupole order (which is sufficiently accurate for 
our purpose), 
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where u = irMf. 
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Fig. 3. — Root noise spectral density for LISA 

In this paper, we shall adopt sensitivity curves 
for LISA developed independently by Larson, His- 
cock & Hellings (2001) and Armstrong, Estabrook 
& Tinto (1999) (sec figure 3). The two methods 
are in substantial agreement, and the former has 
been summarized in the Sensitivity Curve Gen- 
erator (SCG), available online 2 . The sensitivity 
curves incorporate sources of instrumental noise in 
LISA, such as laser shot noise, acceleration noise, 
thermal noise, etc., coupled with a LISA transfer 
function which takes into account the effect of the 
finite time of propagation of the laser beams dur- 
ing the passage of the gravitational waves. We 
assume the case where all three arms of LISA 
are of equal length, and we assume that aver- 
ages over angles and polarizations have been done. 
The SCG permits various choices to be made for 
LISA instrumental parameters, and has an option 
to include an estimate for confusion noise resulting 
from a background of galactic white-dwarf binaries 
(Bender & Hils 1997); this background is included 
in the analysis, but in fact plays no significant role 
for the late stage of IMBH inspiral because the 
signals are predominantly at high frequency, well 
above the white-dwarf band. 

We consider IMBH systems containing one 
black hole in the mass range 10 - 100 M Q and 
a companion black hole in the mass range 100 - 
1000 Mq. Figure 1 shows the results for the dis- 
tance as a function of reduced mass parameter for 
various total masses, assuming a SNR of 10 with a 
one-year integration time ending at the ISCO (or 
at the termination of the LISA sensitivity curve, 
whichever comes first). The distances, in Mpc, 
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can be fit by the approximate formula 




This scaling can be understood as follows: the 
waves from IMBH in the final year of inspiral are 
mainly at the high frequency end of the LISA 
sensitive band, where S n (f) oc f 2 . Also, for 
these masses and year-long integrations, the sec- 
ond term in equation (9) dominates, so that /j = 
(1/8tt)(5/T) 3 / 8 A1- 5 / 8 < //. Consequently' the 
integral B(7) ~ ./T 10/3 ~ ^25/i2 T 5/4 Substi _ 
tuting this into equation (7), we obtain the scaling 
in M. shown in (10). A T 1 / 2 integration-time scal- 
ing actually fits the curves somewhat better than 
a T 5// ® scaling inferred from the analytic estimate. 
Equation (10) underestimates the distances at the 
low mass end, e.g. for 77 <~ 1/10 and M ~ 100 M©, 
by about 20 percent. 

It is straightforward to show that the distance 
reached is independent of LISA acceleration noise, 
but is inversely proportional to the LISA position 
noise error budget. Halving the position noise 
doubles the distance reached. The distance is 
weakly dependent on LISA arm length, varying 
by between 10 and 40 percent for factor-of-two 
changes in arm length; varying arm length does 
not alter the overall level of the high-frequency 
end of the noise curve, but instead moves the lo- 
cation of the peaks in the oscillations of the trans- 
fer function and raises the floor of the noise curve 
near 10~ 3 Hz. 

We obtain similar results using a sensitivity 
curve derived by Finn & Thorne (2000) based on 
work of the LISA Mission Definition team (see ref- 
erence 44 of Finn & Thorne (2000)). This curve 
closely matches the curve from the SCG except for 
the oscillations in the noise root spectral density 
at high frequency resulting from the LISA transfer 
function, an effect of the finite arm length. The os- 
cillations have the effect of reducing somewhat the 
distances inferred from the SCG for a given SNR 



compared to those from the Finn-Thornc curve, 
for sources at high frequency. 

We also calculate the distances that can be 
reached in one year's integration for IMBH sys- 
tems much earlier than the ISCO stage (see Figure 
2). All systems detected within about 40 years of 
merger, can be detected at the Virgo cluster (19 
Mpc) with SNR=10 and one year of integration. 
However at earlier epochs the reach decreases dra- 
matically, so that by 400 years no systems can be 
seen at Virgo, and by 1000 years before merger, 
the distance reached for SNR=10 is in the range 
5.6 to 7 Mpc for all IMBH mass ranges consid- 
ered. This is well short of the Virgo cluster. For 
systems detected 10 6 years before merger, the dis- 
tance reached is 7 Kpc, independent of mass. This 
behavior can be understood analytically. For ob- 
servation epochs long before the ISCO, one can 
show that the inital and final frequencies observed 
are related by w //(l + T/T')~ 3/s , where T 
is the integration time, and T" is the epoch before 
merger, with T <C T 1 . Then, in a region where the 
noise spectral density is dominated by the behav- 
ior S n (f) w a/", one can show, using equations 
(7) and (8) that 

(ii) 

For epochs earlier than 1000 years before merger, 
the systems are in the low-frequency regime where 
S n (f) oc / (apart from white-dwarf confusion 
noise), hence D L w (l/32 7 r 2 )(T/a) 1 / 2 (TV)- 1 , in- 
dependent of chirp mass. With a = 6.09 x 
10~ 51 Hz 3 giving a good fit to the low-frequency 
LISA instrumental noise, we obtain 

D L w 7T 1/2 (1000/T')(10/p)Mpc, (12) 

which perfectly matches the large-epoch behavior 
in Figure 2. The T 1 / 2 behavior is expected for 
what are observations of an essentially CW source. 
These results also are far more pessimistic than 
the estimates made by Miller & Colbert (2004), 
which suggested detecting IMBH binaries in the 
Virgo cluster 1000 years before merger, and in the 
local system of globular clusters 10 6 years before 
merger. On the other hand, because of the T 1 / 2 
behavior, such binaries could be reached in an ex- 
tended LISA mission that enabled 10-year integra- 
tion times. 
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3. Rate of detectable mergers 

To calculate the rate of inspirals detectable by 
LISA in the last year before merger we follow the 
method outlined by Miller (2002). The rate is 
given by 

R=yJ D(Mfu(M)n sc f(M)dM, (13) 

where D{M) is the distance reached for a given 
SNR and integration time, as a function of to- 
tal mass (holding reduced mass fi fixed); v(M) = 
10 _10 ^i _1 M yr _1 is the rate at which smaller black 
holes merge with black holes of mass M in a 
given cluster; n gc = 8h 3 Mpc~ 3 is the number 
density of globular clusters in the local universe; 
and f(M) = [/tot/ln(M max /M min )]M- 1 is the 
fraction of globular clusters harboring black holes 
with mass M per mass interval dM , in the range 
M min < M < M max , with Jf(M)dM = f tot . We 
assume M max » M min . Substituting these formu- 
lae, along with equation (10) into equation (13), 
we obtain equation (2). 

4. Discussion 

Both the distance reached by LISA and the es- 
timated rate of inspirals detected differ markedly 
from the estimates given by Miller (2002), namely 
D L w 2OO(m/1OM ) 1 / 2 (M/1OOM ) 1 /2m P c, and 
equation (1). The distance derived by Miller 
is larger and the mass dependence is different 
than in equation (10) because he appears to have 
adopted an equation [(B7)] from Flanagan & 
Hughes (1998) that actually applies only to equal- 
mass systems. For the low masses relevant to 
IMBH in the year leading up to the merger, the 
high-frequency end of the LISA noise spectrum 
is dominant, where the dependence on frequency 
(and hence on M) is steeply increasing. But the 
underlying dependence is on chirp mass, so for a 
given total mass, there is still a strong dependence 
on reduced mass parameter, which is not reflected 
in equation [(B7)] of Flanagan & Hughes (1998). 
Also, the noise root spectral density curve mod- 
eled by Flanagan & Hughes (1998) is lower by a 
factor of about 3-5 than that given by the SCG 
(figure 3). Their noise curve was based on an out- 
of-date description of the LISA mission. Distance 
estimates using the Flanagan-Hughes curve will 



therefore automatically be larger and rates will be 
higher than those using the SCG. 

We have assumed throughout that the inspi- 
ral orbits are quasi-circular, that is, circular apart 
from the adiabatic inspiral due to radiation reac- 
tion. In reality, IMBH binary orbits are likely to 
be highly eccentric (Gultekin, Miller & Hamilton 
2004). This will increase the average gravitational- 
wave flux for a given orbital period, and could add 
harmonics in frequency bands where LISA's noise 
may be lower; on the other hand it will increase 
the number of parameters to be estimated in the 
matched filter, which raises the threshold needed 
to achieve detection. It is not clear whether the 
net effect of these complicated and possibly offset- 
ting effects will improve LISA's reach for IMBH 
binaries or not. In any event, they are beyond the 
scope of this paper. 
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